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ABSTRACT 

A series of five unusual slow glitches of the radio pulsar B1822-09 (PSR J1825-0935) 
were observed over the 1995-2005 interval. This phenomenon is understood in a solid 
quark star model, where the reasonable parameters for slow glitches are presented in 
the paper, ft is proposed that, because of increasing shear stress as a pulsar spins down, 
a slow glitch may occur, beginning with a collapse of a superficial layer of the quark 
star. This layer of material turns equivalently to viscous fluid at first, the viscosity 
of which helps deplete the energy released from both the accumulated elastic energy 
and the gravitation potential. This performs then a process of slow glitch. Numerical 
calculations show that the observed slow glitches could be reproduced if the effective 
coefhcient of viscosity is ^ 10^ cm^-s~^ and the initial velocity of the superficial layer 
is order of 10~^^ cm-s~^ in the coordinate rotating frame of the star. 

Key words: dense matter - stars: neutron - pulsars: general - pulsars: individual: 
PSR B1822-09 



1 INTRODUCTION 

It is a pity that we are now not able to determine confidently 
which state of matter really exists in pulsar-like stars due to 
the difficulty of calculation from the first principles of the 
elementary strong interaction, even 40 years after the dis- 
covery of pulsars. Nuclear matter (related to neutron stars) 
is one of the speculations even from the Landau's time, while 
quark matter (related to quark stars) is an alternative due to 
the fact of the asymptotic freed om of strong interaction be- 
twccn quarks (see reviews, e.g., 'Madsen"l999': i Glendenningj 
1.2000; Lattimcr & Prakash 2004; Weber 2005). We then have 
to focus on astrophysical observations in order to solve this 
important and fundamental question. 

Actually, pulsars are ideal astro- laboratories for physics 
of cold matter at supranuclear density. Based on the Planck- 
like spectrum without atomic features and the precession 
properties of pulsar- like stars , a solid quark matter conjec- 
ture was suggested l|Xull200^ . Additionally, other features 
naturally explained within this mod el could possibl y include 
sub-pulse drifting in radio emission (|Xu et al.1119991 ). glitches 
l|Zhou et al.ll2004l ). strong magnetic fields, birth after a suc- 
cessful core-collapse supernova, and dete ction of s mall bolo- 
metric radii (for a short review, see, e.g., l5Cul|200 6h. Besides, 
the observational features of anomalous X-ray pulsars/soft 
gamma-ray repeaters may also reflect the nature of solid 
quark stars l|Qwenll2005l : iHorvathllioOsI : IXu et al.ll2006l ). 

The slow g litch es re cently observed in P S R B1822- 
09 JShabanw a 1998; Shab anova fc Uramal [20051 : IZou et all 
l2004l : IShabanovai2005 . .20071 ) could be a new probe into cold 
solid quark matter. The pulsar has a period of 0.769 s and 



a relative young age of 230 kyr. What is interesting in this 
pulsar is the change in its period known as slow glitches. 
Different from the typical feature of glitch phenomenon of 
other pulsars, this pulsar shows glitches that have rather 
slow increase in spin frequency, with a time scale of 200-300 
days (while the normal glitches experience this process much 
swifter, maybe less than a spin period), and moreover, this 
pulsar didn't experience any relaxation progress during the 
post slow glitch days. 

A series of five slow glitches are shown in Fig.l. Fig. 
1(a) shows the frequency derivative !>. The peaks of z> are 
enveloped in a parabolic curve, that may indicate that all the 
slow glitches could be the components of one process. This 
process could be triggered by the small glitch {Av/f ~ 8 x 
10"^") occurred in 1994 September. Both Fig. 1 (b) and (c) 
show the frequency residual Av, but relative to fits to data 
for different intervals (1991-1994 for the former, and 1995- 
2004 for the latter). In this observational result, it is evident 
that the typical character of slow glitch in Fig. 1(b): the spin 
frequency experiences a process that increases gradually but 
never decreases (i.e. no post-glitch relaxation there). 

We are trying to explain this phenomenon in the regime 
of solid quark stars. Since the pulsar is supposed to be mostly 
built up of solid quark matter and the typical density there is 
order of 10^* g/cm"^, soon after a quake, the motion of only a 
thin layer of the quark matter may cause significant change 
in the moment of inertia (and thus the spin frequency) of 
the star owing to the extremely high density. Numerical cal- 
culation is carried out to simulate the observational results 
of both the slow increase in the rotating frequency and the 
lack of the relaxation of the post-glitch progress. 
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Figure 1. The observational data of frequency residual Au and 
frequency derivative u for the pulsar B1822-09. Arrows point- 
ing downwards indicate the epochs at which five slow glitches 
occurred, while arrows pointing upwards indicate two normal 
glitches. The figure is drawn by Shabanova (2007). 



2 THE MODEL 

We suggest that the star is a solid quark star with a typical 
density of lO^^g/cm^. Due to magnetodipole radiation and 
particle ejection of the star, its rotation energy loses and 
frequency of the rotation declines. While on the other hand, 
provided that the pulsar cools down from a hotter fluid el- 
lipse, the equilibrium of the pulsar, with spin frequency i^, 
results in an eccentricity of e. 



2 o 3^ 

y — Stv Gp[ 



- (3 — 2e ) arcsin e 



3(1 -e^), 



(1) 



where G is the gravitational constant and p is the aver - 
aged density of the star (e.g., IShapiro fc Teukolskvlll983 ). 
The configuration of Eq.([l| is for spinning stars made up of 
homogenous fluid matter. Quark stars could be well approx- 
imated by unif orm density if their masses a re not greater 
than ~ 1.5M0 (|Alcock. Farhi fc Olintolll986l '). Eq.([T| could 
be simplifled for stars with small eccentricities as long as 
they spin not very fast. 



1/ = 4:Tve 



2TvpG 
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(2) 



The relation of Eq.(l2]) implies that when the frequency 
of a star decreases, the eccentricity of the star should de- 
crease accordingly in order to maintain the equilibrium 
shape if the star is in a fluid state. But, as the pulsar cools 
down and stays as solid star, the elastic stress might oc- 
cur and accumulate to resist the changing of the star shape. 
When this stress develops to a critical value that the ma- 



terial of the star cannot afford, part of the star breaks (or 
collapses) and then the elastic e nergy releases. This is know n 
as a quake of solid quark star (jZhou et al.ll2004l : I5cull2006l l. 
during which both elastic a nd gravi t ationa l energies are re- 
leased. In the model of .Zhou et al.l l|2004l ). the moment of 
inertia, /, decreases suddenly (to result in a sharp jump 
of Ai/), and then increases gradually (to result in a follow- 
ing post-glitch relaxation). However, as elastic force devel- 
ops to a critical value, a solid rotator with smaller shear 
modulus might not be so violent (i.e., I does not decrease 
suddenly) but in a gentle way for / to decrease, and could 
thus reproduce the observational feature of slow glitches. 
No J-incensement occurs if no signiflcant elastic energy is 
converted to kinematic energy (and thus no post-glitch). 
The smaU glitch at 1994 of PSR B1822-09 may lead to a 
small effective shear modulus (e.g., signiflcant part of the 
quark matter would be shiver-like) , and thus trigger the slow 
glitches. A pre-glitch could then be expected for slow glitches 
in the star-quake model of quark stars, in order to have an 
effective small shear modulus. 

Soon after a quake, we assume that a superflcial layer 
moves in order to set a new equilibrium (see Fig. 2), while 
most of the inner part of star keeps almost the same. Part 
of the energy released turns into heat and melts the debris, 
while the other part turns to the kinetic energy of the fluid. 
It is suggested in this model that the broken material will 
turn into viscous fluid and is therefore able to flow slowly 
and change the shape of the star towards the equilibrium 
shape. This progress would surely cause slow decreasing of 
the moment of inertia of the star. Since this process happens 
in a pretty short time comparing to the typical age of the 
pulsar, it is reasonable to suppose that the angular momen- 
tum of the star keeps invariant during the whole progress of 
glitch. As a result, the frequency of the star might increase 
slowly during the same time. When the fluid is flowing, the 
viscous interaction among the elements of the fluid may ex- 
haust the kinetic energy of the fluid. The shape of the star 
stop changing flnally, and the material then cools down and 
is solidifled. This is why slow glitch appears no relaxation 
feature in the post-glitch progress. Calculational results (see 
§3) explain why this pulsar experiences a serious of flve slow 
glitches rather than a single one. 

To assure that model, we are doing numerical calcula- 
tion to simulate this process. First, the Navier-Stokes equa- 
tion is known as the basic equations to describe the motion 
of viscous fluid. In the frame which is flxed on the inner solid 
part of the star, the general form reads, 



d , 1„ . Mlu 

—V + (v ■ V)v = --Vp + pA-v + g - (— X r + 

ot p at 

+UJ X (cj X r) — 2v X oj). 



(3) 



where v is the velocity of the superflcial matter in the rotat- 
ing frame, p is the pressure in the star, g is the acceleration 
of gravity on the surface of the star, p denotes the viscos- 
ity of the fluid after the phase transformation (from solid 
to fluid-like matter), and the star is assumed to rotate con- 
stantly. The three terms in the bracket of right hand side 
arises from the inertial acceleration due to the non-inertial 
frame we choose. 

Let's consider briefly Eq.(|3]). For a zero order approx- 
imation, a solid star could be modelled by a rigidity body. 
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Figure 2. Sketch of the model for slow glitches. The left side of 
this figure shows the change of the eccentricity of the star during a 
slow glitch and the layer activity in this progress. The right side of 
this figure shows the change of u and / versus time, respectively 
(Note: the usual change of u and I due to the energy lost by 
magnetodipole radiation and particle emission is neglected here). 
A slow glitch is supposed to occur at time ti. In the left, (1) an 
imagined ellipsoidal figure (determined by Eq. JTJ) at ti if the star 
is in a fiuid state, with spin frequency u = i-'bcfore; (2) the real 
figure of a solid star at ti (y = I'before); with stress energy high 
enough to quake; (3) the figure after slow glitch at time ta, with 

V = i^afton which might be close to the imagined shape (1); (4) 
the inner solid part which keeps almost the same during a slow 
glitch; (5) the superficial crust layer which is transmuted during 
the slow glitch progress. 

with velocity v = 0, i.e., 

= — ivp + g — (4- X r + X ('^ X r)). (4) 

p at 

The above equation could at least be adaptable for the case 

01 V <^ ruo. Observationally, the timescale of slow glitches is 
much longer than the spin period. We thus think that Eq.Q 
could be applied in our following simulation. When a glitch 
happens, the elastic energy is released suddenly, the solid su- 
perficial layer turns to a fiuid-like state with strong viscosity, 
but the pressure gradient Vp may keep to be almost invari- 
ant. It is a key point that Vp remains approximately the 
same in spite of stress release in the model, which means 
that this is the case in which most of the elastic energy 
released is changed into heat, rather than into the kinetic 
energy. The heating may melt down solid quark matter to 
become viscous fluid. Soon after the starquake, the matter 
could be re-solidified finally due to cooling. Contrastively, in 
case that the star breaks globally, most of the elastic energy 
may turns into the form of the kinetic energy, and the pres- 
sure gradient would change significantly during starquake; 
the star would experience then a normal glitch with sharp 
frequency increase. Therefore, during a glitch, only uj and 

V in Eq.@ change, with Eq.Q to be well approximated 
since Atj/tj ~ 10~^ and Aiii/ti; ~ 10~^ are very small. Then 
Eq.([3ll of the viscous fluid can be reduced to be 

+ (v • V)v = ^iAv - (-2v X cj). (5) 

We are concerning the initial and boundary conditions 
for our simulations. First we assume that only part of a 
solid quark star breaks under the stellar surface (the shaded 
region in Fig. 2), with a thickness of H . We note that this as- 
sumption is for simplifying the calculation, but is not a key 
for general results in the starquake model for slow glitches 
(see §3). We denote h to measure the height of matter in 



the beaked part, i.e., /i = in the bottom while h — H on 
the surface. Two cases are considered in the simulated re- 
suh of Fig. 3, H = 10%R and 20%R, where R is the steUar 
radius. We suppose that the velocity of matter increases lin- 
early from the bottom of collapsed layer to the top, with a 
boundary of zero velocity at the bottom, and that the veloc- 
ity increases gradually from polar to equator. The tangential 
velocity could then be formulated as 

J «max(6') = ViSinO, 

\ v{e,h) = wW-l, ^ ' 

where Vi is the maximum of the initial velocity soon after 
a quake. A starquake is assumed to be a trigger for the 
initial movement with maximum velocity Vi. The continues 
condition of matter (i.e. V • v = 0) determines the radial 
velocity, given that the viscous fluid is incompressible. 

Solving Eq.lO with boundary Eq.(|S}, one can obtain 
a 3-dimensional vector field of velocity all over the star. To 
achieve the changing of the shape and rotating frequency, we 
adopt the incompressible hypothesis and the conservation 
law of angular momentum. The movement of incompress- 
ible matter in the layer changes gradually the shape of star, 
and thus the angular frequency because of conservation of 
angular momentum. In fact the conservation law of angu- 
lar momentum is not obeyed precisely due to the spindown 
torque. However, in hght of the rather short duration of the 
glitch progress, the conservation could be well approximated 
during a relatively short time. 

After a series of slow glitches, we may suppose that 
the star stays at an equilibrium state which is described by 
Eq.([2}. It is worth clarifying that, as noted in Fig. 2 (note 
the difference between states (1) and (3)), after each slow 
glitch (except the last one), the star did not turn to the real 
equilibrium state, but suffer actually a series of glitches, not 
only a single one. We can then obtain the eccentricity of the 
star with the angular frequency at the end of the five slow 
glitches through Eq.((2]). After doing that, we may calculate 
the star's spin frequency in any time via the conservation 
law of angular momentum. 

Getting all above, we can carry out numerical calculus 
after setting the value of initial velocity Vi and viscosity fi. 
Our goal is to find whether there is reasonable parameter 
space (of Vi and fi) where the general features of the slow 
glitches observed could be reproduced. In the calculation, 
the "heun" method is adopted in order to obtain more ac- 
curate values from numerical calculations. 



3 THE RESULTS 

We do the simulation with typical density of p = 10^'' g/cm^ 
and stellar radius of _R = 10 km for indications. However, 
according to various simulations, the general conclusions 
would not change significantly if these two parameters are 
in the same orders. Observational data of slow glitches and 
the simulated result in the solid quark star model are shown 
in Fig. 1(b) and Fig. 3, respectively. The parameters of ^ 
and Vi for the five slow glitches are also shown in table 1. 

In Fig. 3, two simulated curves are shown with differ- 
ent thickness of the surface fiuid layer, H — 10%R and 
H — 20%R, respectively. From Fig. 3, it is obvious that the 
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Figure 3. The simulation of slow glitches. The vertical axis is the 
frequency change, Au (in unit of 10~^ Hz), and the horizontal axis 
measures the time (in days) . The two curvatures in this figure are 
two simulations, given that the thickness of the surface layer to be 
10%R and 20%i?, respectively. The starting point of time could 
be at MJD 49857 in Fig. 1. 
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Figure 4. Different curves of At/ with different parameters of Vi. 
The curve labelled as "Vi" is that same one labelled as "10%R" 
in Fig. 3, with i)i-values listed in Table 1. The curves labelled as 
"2Vi" and "0.5Vi" are with ti^-values double and half the ones 
listed in Table 1, respectively. 
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10-10.88 
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10-10.17 
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102.6 


10-11.10 





Table 1 . The parameters in the simulation where the thickness of 
the surface layer is set 10%_R. The corresponding simulated result 
is shown in Fig. 3. 
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thickness indeed affect the resuh of Ai/. Though the gen- 
eral behavior keeps, the amplitude and the gradient of Av 
increase as percentage of the thickness to the stellar radius 
increases. 

What if the values Vi of and /i in Table 1 change? In 
case oi H = 1Q%R, the influence of these parameters on 
the simulated result of Av are shown in Fig. 4 and Fig. 5, 
respectively. 

As shown in Fig. 4, where values of ^ are the same 
listed in Table 1, Aiz-curves for various Vi are illustrated. It 
is evident that a larger value of Vi could result in a larger 
glitch amplitude of Av. It is worth noting, however, that 
the time durations in which Av increases to the maximum 
are almost the same. This could be easily understood since 
a larger Vi implies a larger stress and energy released in the 
starquake. The effect of changing /.i-values on the amplitude 
of Av are demonstrated in Fig. 5, where the values of Vi are 
the same listed in Table 1. As shown in Fig. 5, smaller fi- 
values could not only result in higher amplitudes of Av, but 
also in longer durations for Av to increase to maximum. This 
result is also understandable. In the model, the fluid moves 
under coriolis' force and viscous resistance. Smaller p could 
lead to a relative smaller resistance and thus a longer time 
to reach its equilibrium. In summary, one needs to have a 
high Vi or a low n in order to reproduce a large amplitude of 



Figure 5. Different curves of Av with different parameters of fi. 
The curve labelled as "/i" is that same one labelled as "10%R" 
in Fig. 3, with /^-values listed in Table 1. The curves labelled as 
"2/i" and "0.5/^" are with /x- values double and half the ones listed 
in Table 1, respectively. 

Av, while a small i^- value could effectively delay the increase 
of Av. 



4 CONCLUSIONS AND DISCUSSIONS 

A solid quark star model, in which the surface matter breaks 
during an initial small glitch, is suggested for understand- 
ing the general features of slow glitches observed recently. 
Though this is a real problem of rheology in principle, we 
deal with the matter as viscous fluid for simplicity. Numer- 
ical calculations show that the observed slow glitches could 
be reproduced if the effective coefficient of viscosity is ~ 10^ 
cm^.s"^ and the initial velocity of the superficial layer is or- 
der of 10~^° cm-s"^ in the coordinate rotating frame of the 
star. We simulated with typical density of p = lO^^ g/cm^ 
and stellar radius of _R = 10 km for indications, but the gen- 
eral conclusions would not change significantly if these two 
parameters are in the same orders. 
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The pieces of blocks should be re-correlated through 
solidifying melted surfaces of the pieces, into which the su- 
perficial layer of a solid quark star breaks soon after a quake. 
The released energy (both gravitational and elastic) during 
a slow glitch may result in melting the conjunctional parts 
of segments (i.e., the surfaces of blocks) at first. The tem- 
perature then increases significantly at only a small fraction 
(77 <^ 1) of the segments. For an order-of-magnitude demon- 
stration, we apply the Debye model to normal solid quark 
matter and obtain a heat capacity of 



Cv = 0.15- 



PB 



-V- 



(7) 



The matter absorbs an energy of E = CvdT from an 
initial temperature of To ^ T to T. In case of released 
energy of i? ~ 10^^ erg and a melting temperature of T ^ 10 
MeV (^ the star's global temperature To ~ keV), we have. 



2.7 X 10" 



10" 



RlT,\ PO ' 



(8) 



for stars with radius of Ro x 10® cm and density of p, where 
Til X 10^^ K is the melting temperature and po is the nuclear 
density. For blocks with length scale of 1 m, this calculation 
shows a heating surface part with thickness of about lO'^ fm 
for each of the blocks. We note the timescale for keeping 
such a high temperature gradient should be very small, so 
that the contacting parts cool rapidly to solidify. More and 
more small blocks become correlated (i.e., a big solid-like 
bulk matter forms) as a slow glitch evolves, and strain energy 
develops again then. All the heated points are within the star 
and the cold surface keeps, and then no significant heating 
feature could be observed after a glitch. 

Alternatively, let's consider a crusted strange quark star 
with a solid crust and a fluid quark matter core. The iner- 
tia of the crust could be only ~ 10~^ times of that of the 
quark matter core. According to AI/I ~ AQ/ft ~ 10~* for 
slow glitches, the crust ellipticity should change a value of 
~ 10~^, while the actual ellipticity is ~ 10~^ for stars with 
period of ~ 1 s. This means that pulsar slow glitches could 
not be re-produced in a crusted strange star model. 

What is the key ingredient which affects a pulsar to 
undergo a normal or a slow glitch? We think that the stellar 
mass of quark star could play an important role in this issue. 

Let's introduce t wo k i nds of stress force inside solid stars 
at first. As noted bv IXul ()2006l ). two kinds of factors could 
result in the development of stress energy in a solid star, 
and then in star quakes as glitches, (i) As a quark star cools 
(even spinning constantly), changing state of matter may 
cause a development of anisotropic pressure distributed in- 
side a solid matter. Such matter cannot be well approxi- 
mated by perfect fluid, and the equation governing star's 
gravitational equilibrium should then not be the TOV equa- 
tion. In case of spherical symmetry (the simplest case), one 
can introduce the difference between radial and tangential 
pressures, A. Change of A would lead to no-conservation of 
stellar volume. We call the force, which acts primary in this 
case, as bulk-variable force. This kind of force may be the 
key factor for giant quak es during superflares of soft 7-ray 
repeaters (|Xu et al.ll2006l) . (ii) An uniform fluid star would 
keep its eccentricity presented in Eq.((T)) or Eq.lQ, i.e., the 
eccentricity decreases as a star spins down. However, for a 
solid star, the shear stress would prevent the star from de- 



creasing eccentricity during spindown. In this case, even the 
state of matter does not change, stress energy could still 
develop as solid star spins down. We call this kind of force 
as bulk-mvariable force since the total stellar volume may 
always keep constantly. Starquakes due to bulk-invariable 
force, and the consequent g litches, were calculated and dis- 
cussed in (|Zhou et al.ll2004l ). 

Both bulk-invariable and bulk- variable forces could re- 
sult in decreases of moment of inertia, and therefore in pul- 
sar glitches. These two kinds of forces could trigger nor- 
mal glitches if they are relatively stronger than the critical 
stress, but might only conduce to slow glitches if weaker. 
The bulk-variable force could be stronger than the bulk- 
invariable one for massive solid stars, since, in a special case 
of non-rotating, the gravity oc K'P /R^ [M: mass, R: radius) 
is strong there. If the critical stress of solid quark matter 
is almost the same (note: the effective critical stress could 
be much small for a shiver- like surface layer), we may think 
that normal glitches prefer to occur in massive quark stars 
with possibly strong bulk-variable forces, while slow glitches 
are likely to take place in low-mass solid quark stars. 

What if PSR B1822-09 is a low-mass quark star? The 
spindown due to magnetodipole torque for a star with a 
magnetic dipole moment p, a moment of inertia /, and an 
angular velocity SI reads 



3Jc3 



2o3 
p it . 



(9) 



This rule keeps quantitatively for any oblique rotators, as 
long as the braking torques due to magn etodipole radiatio n 
and the unipolar generator are combined l|Xu fc Qiacill200il ). 
Approximating I = {2/5)MR^, p = {1/2)BR^, and M = 
{4:7V /S)R'^ p, with p the average density, we have then 



R 



2.76 X lO^VwB"^ cm. 



(10) 



where p = pi4 x lO^^g/cn^, and P = 0.769 s and P = 
5.23 X 10'^" for PSR B1822-09. The pulsar could be a few 
kilometers in radius if the polar magnetic field is ~ lO'^^ G. 

The potential drop in the open field line region, 
should be greater than a critical value of ~ 10^^ V ~ 3 x 10^ 
e.s.u. in order to create secondary plasma via gap dis- 
charges. The potential dro p between the center and the edg e 
of a polar cap is therefore l|Ruderman fc Sutherlandlll975l ) 



^R^BP'^ 
c^ 



7.84 X IQ'^pl^B- 



(11) 



Note that the potential in above equation would be higher if 
the effect of inclination angle is included (|Yue et al.ll2006l ). 
The drop could be high enough for pair production if the 
field B ~ 10^^ G, and the pulsar should then be radio loud. 
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